Unconventional Integer Quantum Hall effect in graphene 
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Monolayer graphite films, or graphene, have quasiparticle excitations that can be described by 

2 + 1 dimensional Dirac theory. We demonstrate that this produces an unconventional form of the 
quantized Hall conductivity a xy — — (2e 2 //i)(2n + 1) with n = 0, 1, . . ., that notably distinguishes 
graphene from other materials where the integer quantum Hall effect was observed. This unconven- 
tional quantization is caused by the quantum anomaly of the n = Landau level and was discovered 
in recent experiments on ultrathin graphite films. 
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The quantum Hall effect (QHE) is one of the most re- 
markable phenomena in condensed matter discovered in 
the second half of the 20th century. The basic experimen- 
tal fact characterizing QHE is that the diagonal electric 
conductivity of a two-dimensional electron system in a 
strong magnetic field is vanishingly small a xx —> 0, while 
the non-diagonal conductivity is quantized in multiples 



of e 2 /h: 
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-ve /h, where v is an integer (the integer 



quantum Hall effect (IQHE)) or a fractional number (the 
fractional QHE). In a recent paperi the fabrication of 
free-standing monocrystalline graphite films with thick- 
ness down to a single atomic layer was reported. This 
new material, called graphene, possesses truly remark- 
able properties such as excellent mechanical characteris- 
tics, scalability to the nanometer sizes, and the ability 
to sustain huge (> 10 8 Aj cm 2 ) electric currents. By us- 
ing the electric field effect, it is possible to change the 
carrier concentration in samples by tens times and even 
to change the carrier type from electron to hole when 
the sign of applied gate voltage is altered. All this make 
graphene a promising candidate for applications in future 
micro- and nanoelectronics. 

On the theoretical side, the linear, Dirac-like, spec- 
trum of quasiparticle excitations (up to energies of the 
order of 1000 K) and the pseudospin degeneracy make 
graphene a unique truly two-dimensional " relativistic" 
electronic system. The thinnest graphite films can be 
described by a low-energy (2+1) dimensional effective 
massless Dirac theory^. Of special interest are the prop- 
erties of graphene in a magnetic field. The important 
differences between the Dirac and Schrodinger theories 
may be observed in thermodynamic and magnetotrans- 
port measurementsi^i^i&. For instance, the phase of 
de Haas van Alphen and Shubnikov de Haas oscillations 
for Dirac quasiparticles is shiftec£i2i2ii& compared to the 
phase of non-relativistic quasiparticles. Moreover, the 
Dingle and temperature factors in the amplitude of os- 
cillations explicitly depend on the carrier density in the 
case of a Dirac-like spectrum^. 

Because of the large value of the cyclotron gap, it is 
expected that the QHE in this material can be observed 
for much higher temperatures and lower magnetic fields 
than in conventional semiconductors. Therefore it is nat- 



urally to ask whether the fundamental difference between 
the properties of Landau levels (LL) (see Eqs. © and 
HI U| l below) in the Dirac and Schrodinger theories can be 
observed experimentally in the Hall conductivity? The 
purpose of this letter is to show that the Dirac-like dy- 
namics of graphene results in an unconventional form of 
the Hall quantization 
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We argue that the quantization rule Q is caused by the 
quantum anomaly of the n — LL , i.e. by the fact 
that it has a twice smaller degeneracy than the levels 
with n > and its energy does not depend on the mag- 
netic fielcUi. Remarkably this quantization is observed 
experimentally^ for ultrathin graphite films which ex- 
hibit the behavior expected for ideal 2D graphene. 

We begin with the Lagrangian density of noninteract- 
ing quasiparticles in a single graphene sheet that in the 
continuum limit reads 2 
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where ty a = (ipi a (t,r),ip 2a (t,r)) is the four-component 
Dirac spinor combined from two spinors tpi a , 4>2a [corre- 
sponding to K and K' points of the Fermi surface, re- 
spectively] that describe the Bloch states residing on the 
two different sublattices of the biparticle hexagonal lat- 
tice of the graphene sheet, and a = ±1 is the spin. In 
Eq. 7 M with \i = 0, 1, 2 are 4x47 matrices belonging 
to a reducible representation in 2 + 1, ^> a = is the 

Dirac conjugated spinor, — e < is the electron charge, 
vp is the Fermi velocity. We set ks = 1, but kept Planck 
constant h — h/2n. 

The external magnetic field B is applied perpendic- 
ularly to the plane along the positive z axis and the 
vector potential is taken in the symmetric gauge A = 
(—B/2y,B/2x). In contrast to the truly relativistic 
(3 + 1) casein, the Zeeman interaction term still has to be 
explicitly added to the Lagrangian J5J), because it orig- 
inates from nonrelativistic many-body theory. This can 
be done by considering spin splitting \i a = fi — o-[AbB 
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of the chemical potential p, where ps — eh/(2mc) is 
the Bohr magneton. However, for realistic values of 
vp oc 10 5 m/s in graphene the distance between LL is 
very large compared to the Zeeman splitting^, so that 
in what follows we will not consider this term and sim- 
ply multiply all relevant expressions by 2 to count the 
spin degeneracy. While simple tight-binding calculations 
made for the hexagonal lattice of a single graphene sheet 
predict that p = 0, the real picture is more complicated 
and the actual value of p can be nonzero due to finite dop- 
ing and/or disorder. Moreover, nonzero and even tunable 
value of p [including the change of the character of car- 
riers, either electrons or holes] is possible in electric-field 
doping cxpcrimentsA*^*i£. In our notations p > corre- 
sponds to electrons and accordingly to the positive gate 
voltage. 

Using the Kubo formalism and modeling the LL by 
Lorentzians with a constant width T the following ex- 
pression for the diagonal conductivity was obtained in 
Refs&H 

2e 2 f°° , 
a xx (B,p,,T) = — / <kj[-n' F (w - p)]A xx (uJ, B,T), 

(3) 

where n F (u>) = l/[exp(w/T) + 1] is the Fermi distribu- 
tion and the function A xx that incorporates the effect of 
all LL is given by Eq. (11) of Ref£. Now this result is ex- 
tended for the Hall conductivity and we derive a general 
analytical expression for a xy (B, p, T)m& The resulting de- 
pendence a xy (p) is shown in Fig.^ where one sees that 
the plateaux of a xy follow Eq. • This agrees with the 
latest experimental results^ and resemble earlier theo- 
retical predictions 15 . However, to demonstrate result 
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FIG. 1: (Color online) The Hall conductivity a xy measured 
in e 2 /h units as a function of chemical potential p for two 
different values of F for T = 3K and B = IT. We use 
hv 2 F eB/c -> (4.5 x 10 4 K 2 )B(T). 

in the most transparent way it is useful to write down a 
simpler conventional representation^i^ for a xy obtained 
in the clean limit r — > 0: 
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Here we introduced the filling factor of LL, vb = 
nhc\p\/\eB\ with p being the carrier imbalance (p = 
n e — nh, where n e and rih are the densities of "electrons" 
and "holes", respectively). This filling factor can be rep- 
resented as a sum over LL 
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of the Dirac theory 
1 
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where we separated out the level with n = because its 
degeneracy is only half of the degeneracy of the levels 
with n > 0. To illustrate this rather peculiar property 
of the Dirac theory in a perspicuous way, we included 
in M n and vb the mass (excitonic gap) A which was 
discussed recently to explain some experiments^*^. Our 
consideration of a xy is in fact independent of the presence 
of A, so in what follows we set A = 0. A zero value of 
A is expected for noninteracting quasiparticlcs on the 
hexagonal lattice of graphene. 

The first equality in Eq. (0} corresponds to a clas- 
sical straight line <j\i oc vb- As discussed, for exam- 
ple, in Ref^, this line emerges from two step func- 
tion dependences, viz. p(n) and 012 (/i). Indeed, using 
tanh(w/2T) = sgn(w) for T — > 0, we obtain from Eqs. Q 
and |JBJ (compare with Ref>i&) 
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where [x] denotes the integer part of x. The usual ar- 
gumentation (see e.g. Refii&) for the occurrence of the 
IQHE states that in the presence of disorder the depen- 
dence of (712(1^) remains the same, while p(n) becomes a 
smooth function. The classical (@J and quantum (JJJ) Hall 
conductivities coincide only for the fillings, vb — 2n + 1 
(see Fig. 01 . The odd integer rather than integer fillings 
that produces the quantization rule appears due to 
the above-mentioned halved degeneracy of the n = LL. 
Another interesting feature of Eq. J7|) (see also Figs. 
and|2J) is that a xy — ±2e 2 //i for the fillings vb < 1 and 
it crosses only when p, changes sign. On the contrary, 
in a conventional IQHE a xy — for v B < 1. 

Although Eqs. 0-0 are obtained in the clean limit 
and using a simple bare bubble expression for conductiv- 
ity, our main result is model independent and is only 
based on the n — level anomaly. 

Now we rewrite Eqs. 0-0 in terms of the Fermi dis- 
tribution 
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FIG. 2: (Color online) The Hall conductivity a xy measured 
in e 2 /h units as a function of the filling vb- The straight line 
corresponds to a classical dependence (0J. The parameters 
are the same as in Fig. ^ 



to compare it with Eq. (18) of Ref4£ that was obtained 
for an ideal two-dimensional electron gas 
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with nonrelativistic spectrum 
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There is a commonsense reasoning 1 * 2 * that graphcne 
is a two-band [the first band would corresponds to the 
electrons with to n = M n — [i and the second band, to the 
holes with uj n = —M n — /i], two- valley [corresponding to 
K and K' points of graphen's Fermi surface] semiconduc- 
tor with zero gap A between the bands. Accordingly its 
Hall conductivity can be directly obtained from by 
summing over all these bands and valleys 
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x [nj?(M n — A*) + "f(- M n - fi) 

-n F {M n+ i — fj,) - n F {-M n+ i - fi)] , 
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where we also counted spin degeneracy. It is easy to see 
that Eqs. JSJ and ljTT)l correspond to two completely dif- 
ferent Hall conductivity quantization rules, viz. Eq. (JHJ 
which correctly counts the degeneracy of the n = level 
produces Eq. (JTJ, while the semiconducting analogy l|ll|) 
leads to 
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Here we assumed that e,-B,/i > 0. Although previous 
experimental observations supported the picture based 
on Eq. (|12|) . the latest experiments made on thin films 10 



are in accord with the unconventional Hall quantization 
0. This shows that in an applied magnetic field the 
semiconducting interpretation of graphene's band struc- 
ture that led us to Eq. Ijllfl becomes invalid (see also 
RefiSi). The drastic difference between Eqs. JIJ and i(T2"|) 
is caused by the above-mentioned fact that the lowest LL 
in Dirac theory is special and has twice smaller degen- 
eracy than the levels with n > 0, because depending on 
the sign eB it is occupied either by electrons or holes, 
while higher levels contain both electrons and holes** 2 *. 
In the nonrelativistic theory when the Lande factor g ^= 2 
all Landau levels have the same degeneracjii*. It turns 
out that graphene for which the valence and conduction 
bands intersect in discrete points*^, is reasonably well de- 
scribed by the Dirac formalism which naturally embodies 
the n = level anomaly. 

We now consider the phenomenon of quantum mag- 
netic oscillations in graphene which is closely related to 
the quantization of a xy and discuss the specific of the 
n = level. The de Haas van Alphen and Shubnikov 
de Haas effects in graphene were studied in Refsii**. In 
particular, in Ref£ it was shown that the oscillatory part 
of the diagonal conductivity |j2J) is given by 
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R T (k)R D (k)R s (k), (13) 



where Rt, Rd and R s are respectively the temperature, 
Dingle and spin factors. Using the relationship [i 2 = 
■jrh 2 v F \p\ valid for T = T = B = one can check 
that the minima of the diagonal conductivity occur 
at the fillings Vb = 2n + 1 giving an indication of the 
possible positions of the plateaux in the IQHBiS. [Note 
that in thick films the minima of a xx occur at integer 
fillings 4 .] Obviously for fi = there is no oscillations of 
<j xx , the conductivity <J xx (p = 0) = 2e 2 /(7r 2 /i) becomes 
a field independent universal 2 * quantity that is another 
distinctive feature of the n — level anomaly. 

Although the quantization can be understood by 
considering nonintcracting Dirac quasiparticles placed in 
an external magnetic field, even this simple model reveals 
other unusual properties 11 intimately related to nontriv- 
ial dynamics of quasiparticles from the n = level. For 
example, the t/(4) symmetry of the Lagrangian (J2J is 
spontaneously broken down to U(2) x £7(2) at fj, = 
in non-zero magnetic field even in the absence of ad- 
ditional interaction between fermionsii thus leading to 
the emergence of the chiral condensate (^^). Includ- 
ing many body effects such as an attractive interaction 
between quasiparticles, could further generate a gap for 
quasiparticles like the above mentioned gap A (see e.g. 
Refsfiiftii). Fortunately in the case of the IQHE the pres- 
ence of the condensate does not affect our consideration. 
On the other hand, a possible gap generation for the 
fcrmions from the lowest LL might become important for 
the fractional quantum Hall effect and this issue certainly 
deserves further experimental and theoretical study. 

To conclude, we have shown that the integer numbers 
associated with quantized Hall conductivity in graphene 
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have an unusual pattern a xy h/e 2 = 2, 6, 10, 14 . . .. We 
argued that it is related to the fact that a theoretical 
description of graphene is based on 2 + 1 dimensional 
Dirac theory, where the lowest Landau level has half of 
the higher Landau levels degeneracy. 

We are indebted to A. Geim for showing us his ex- 
perimental results prior to publication and for stim- 



ulating discussions. We also thank J. P. Carbotte, 
V.M. Loktev and V.A. Miransky for useful discussions 
and W.A. de Heer, P. Kim for informing us about their 
latest results. S.G.Sh. was supported by the Natural 
Science and Engineering Council of Canada (NSERC) 
and by the Canadian Institute for Advanced Research 
(CHAR). 



* Electronic address: vgusynin@bitp.kiev.ua 
' Electronic address: sharapov@bitp.kiev.ua 

1 K.S. Novoselov, A.K. Geim, S.V. Morozov, D. Jiang, 
Y. Zhang, S.V. Dubonos, I.V. Grigorieva, and A. A. Firsov, 
Science 306, 666 (2004); K.S. Novoselov, A.K. Geim, 
S.V. Morozov, S.V. Dubonos, Y. Zhang, and D. Jiang, 
cond-mat/0410631 

2 G.W. Semenoff, Phys. Rev. Lett. 53, 2449 (1984). 

3 Y. Kopelevich, J.H.S. Torres, R.R. da Silva, F. Mrowka, 
H. Kempa, and P. Esquinazi, Phys. Rev. Lett. 90, 156402 
(2003). 

4 S.V. Morozov, K.S. Novo selov, D. Jiang, A.A. Firsov, 
S.V. Dubonos, A.K. Geim, cond-mat/0505319 

5 C. Berger, Z. Song, T. Li, X. Li, AY. Ogbazghi, R. Feng, 
Z. Dai, A.N. Marchenkov, E.H. Conrad, P.N. First, and 
W.A. de Heer, J. Phys. Chem. B 108, 19912 (2004). 

6 Y. Zhang, J. P. Small, M.E.S. Amori, and P. Kim, 
Phys. Rev. Lett. 94, 176803 (2005). 

7 S.G. Sharapov, V.P. Gusynin, and H. Beck, Phys. Rev. B 
69, 075104 (2004). 

8 V.P. Gusynin and S.G. Sharapov, Phys. Rev. B 71, 125124 
(2005). 

9 LA. Luk'yanchuk and Y. Kopelevich, Phys. Rev. Lett. 93, 
166402 (2004). 

A.K. Geim, private communication. 

1 V.P. Gusynin, V.A. Miransky, and LA. Shovkovy, 
Phys. Rev. Lett. 73, 3499 (1994); Phys. Rev. D 52, 4718 
(1995). 

2 Note that in the relativistic (3 + 1) Dirac theory the 
Zeeman term is built in the formalism. When the rel- 
ativistic eigenenergy, E(n,o — ±, k z — 0) = [m 2 c 4 + 
mc 2 h\eB\/ (mc)(2n + 14- o")] 1,/2 is separated into the Pauli 
term with afisB and the diamagnetic term, the nonrel- 
ativistic LL with the same degeneracy and the spectrum 
given by Eq. 1101 emerge. Each nonrelativisitic LL suffers 
a Zeeman splitting 4zl/2g/^s_B, so that it doubles the num- 
ber of LL. Here we included the Lande factor g which plays 
an important role in solids. For g = 2, there is an acciden- 
tal (from a nonrelativistic point of view) degeneracy of LL, 
so that each level with n > is filled by both spin down 
and spin up electrons from the n — 1 level. The n = LL 
is special and remains nondegenerate, so that we come to 
the picture considered in this paper. In our case, however, 
an "accidental" degeneracy of n > LL is not related to 
the spin degree of freedom and occurs due to the presence 



of two sublattices in graphene. 

13 E.V. Gorbar, V.P. Gusynin, V.A. Miransky, and 
LA. Shovkovy, Phys. Rev. B 66, 045108 (2002). 

14 The corresponding analytical expression for A xv (lo, B, T) is 
simple, but rather lengthy, V.P. Gusynin and S.G. Shara- 
pov, in preparation. 

1 Similar figures also appear in Y. Zheng and T. Ando, Phys. 
Rev. B 65, 245420 (2002), but with a twice smaller size of 
the steps of o xy , and without making a link between the 
unusual behavior of a xy and the quantum anomaly of the 
n = level. 

16 A.M.J. Schakel, Phys. Rev. D 43, 1428 (1991). 

17 D.V. Khveshchenko, Phys. Rev. Lett. 87, 206401; ibid. 87, 
246802 (2001). 

18 M. Jan/Jen, O. Veihweger, U. Fastenrath, and J. Hajdu, 
Introduction to the Theory of the Integer Quantum Hall 
Effect, edited by J. Hajdu, VCH, Weinheim, 1994. 

19 M. Jonson and S. M. Girvin, Phys. Rev. B 29, 1939 (1984). 

20 M.S. Dresselhaus and G. Dresselhaus, Adv. Phys. 51, 1 
(2002). 

21 T. Ando, J. Phys. Soc. Jpn. 74, 777 (2005). 

22 M.H. Johnson and BA. Lippmann, Phys. Rev. 76, 828 
(1949). 

23 The points near which the electrons have a linear spectrum 
and are described by a two-level Hamiltonian are called 
in the literature— diabolic points . There is a theorem 
of H.B. Nielsen and M. Ninomiya, Nucl. Phys. B185, 20 
(1981); B 193, 173 (1981) that asserts that in condensed 
matter systems diabolic points come in parity-inavariant 
pairs. Accordingly the Dirac Lagrangian J5J which embe- 
des a pair of such points, written using a parity preserving 
4x4 reducible representation of 7 matrices is rather nat- 
ural. 

24 This universality has the same origin as the universality 
considered in E. Fradkin, Phys. Rev. B 33, 3263 (1986) 
for degenerate semiconductors and in PA. Lee, Phys. 
Rev. Lett. 71, 1887 (1993) for d-wave superconductors. 
Moreover in the case of graphene where the quasiparticles 
have not only linear disperion relation similar to the Bo- 
golyubov quasiparticles in d-wave superconductor, but also 
couple to the vector potential in the conventional QED2+1 
way, the conductivity o xx (/j, = 0) is expected to be field 
independent. 



